Towards the electron EDM search: Theoretical study of HfF 
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We report first ab initio relativistic correlation calculations of potential curves for ten low-lying 
electronic states, effective electric field on the electron and hyperfine constants for the 3 Ai state of 
cation of a heavy transition metal fluoride, HfF + , that is suggested to be used as the working state 
in experiments to search for the electric dipole moment of the electron. It is shown that HfF + has 
deeply bound 1 E + ground state, its dissociation energy is D e — 6.4 eV. The 3 Ai state is obtained 
to be the relatively long-lived first excited state lying about 0.2 eV higher. The calculated effective 
electric field E e g = Wd|fi| acting on an electron in this state is 5.84 x 10 24 Hz/(e-cm). 



Introduction. The search for the electric dipole mo- 
ment (EDM) of the electron, d e or eEDM below, remains 
one of the most fundamental problems in physics. Up 
to now only upper limits were obtained for \d e \. The 
tightest bound on d e was obtained in the experiment on 
the Tl atom which established an upper bound of 
\d e \ < 1.6 x 10~ 27 e-cm (e is the charge of the electron). 
It seems unlikely that experiments on atoms can give 
tighter bounds. Molecular systems, on the other hand, 
provide a way for much enhanced sensitivity, since the ef- 
fective intramolecular electric field acting on electrons in 
polar molecules can be five or more orders of magnitude 
higher than the maximal field available in the laboratory 

The new generation of the eEDM experiments, employ- 
ing polar heavy-atom molecules, is expected to reach sen- 
sitivity of 10" 30 -10" 28 e-cm/Vday (e.g., see @|). Their 
results are expected to dramatically influence all the pop- 
ular extensions of the standard model, in particular su- 
persymmetry, even if bounds on the P, T-odd effects com- 
patible with zero are obtained (see [g, 0] and references 
therein). These studies include the beam experiments 
carried out on YbF molecular radicals by E. Hinds group 
at Imperial College, London [5|, and the vapor cell exper- 
iment on the metastable o(l) state of PbO prepared by 
the group of DeMille at Yale University (see [1] and ref- 
erences therein). New ways of searching for eEDM, using 
trapped cold molecular cations, were investigated during 
last two years at JILA by E. Cornell group. The first 
candidate was HI + [§], but subsequent estimate [1(3] and 
accurate calculations have showed that the cation 
has rather small E c g (see below) . Then some more moti- 
vated candidates were considered, HfH + [ijj], HfF + [l3| 
etc., having, possibly, working 3 Ai as the ground or at 
least long-lived first excited state. Even experimental 
study of spectroscopic properties for this cation is a very 
difficult problem and there are no such data measured up 
to date though they are required already to analyze basic 
stages of the eEDM experiment [l3| . In turn, the modern 
relativistic computational methods can now give reliable 
answers on almost all the questions of interest even for 
compounds of heavy transition metals as is in our case. 
So, the goal of the present paper is a theoretical study 



of required spectroscopic properties, E e g and hyperfine 
structure of HfF + as a candidate for the eEDM experi- 
ment. Such a complete analysis is the first one performed 
for the given kind of molecules. 

E e f[, hyperfine structure and working state in HfF^. 
One of the most important features of such experiments 
is that the knowledge of the effective electric field, E c s , 
seen by an unpaired electron is required for extracting 
d e from the measurements. E e g can not be obtained in 
an experiment, rather, electronic structure calculations 
are required for its evaluation. It is presented by the 
expression E e g = Wd|fi|, where Wd is a parameter of the 
P,T-odd molecular Hamiltonian that is given in Refs. [1, 



W d 
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where VP is the wavefunction for the working 3 Ai state, 
and fl = (^{J • n\^f) = ±1, J is the total electronic 
momentum, n is the unit vector along the molecular axis 
directed from Hf to F, 
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E is the inner molecular electric field, er are the Pauli ma- 
trices. Recently, the E e s value for the 3 Ai state was esti- 
mated in the scalar-relativistic approximation by Meyer 
et al. [16( and the method used for calculation of E e s is 
close in essence to that developed by us earlier [TtJ and 
applied to first two-step calculations of the PbF molecule 
1J, [l7( . In the present work, a more reliable value of E e s 
is calculated using the advanced two-step techniques de- 
veloped by us later 18, TjJ 2^, [2l|. The hyperfine con- 



stants for 177 Hf and 19 F nuclei (Ay [Hf] and A|,[F]) for 
the 3 A! state of HfF + are also calculated. When the 
experimental value Ay [Hf] is measured, it will provide 
an accuracy check for the calculated E e g value. Both 
Ay [Hf] and A\\ [F] values arc useful for identifying HfF + 
by his spectrum from another species in the experiment. 
Moreover, preparation of the HfF + cations in the work- 
ing state for the eEDM experiment, registration of eEDM 



2 



signals etc. requires the knowledge of spectroscopic prop- 
erties. In particular, an open question is to clarify, which 
state is the ground one. If 3 A! is not the ground state, 
its radiative lifetime is also required. In the present work 
these data are obtained, though precise studying of the 
compounds with transition elements is a difficult problem 
for modern molecular theory. The needed spectroscopic 
information cannot be obtained now from other sources. 

Methods and calculations. The 12-electron general- 
ized relativistic effective core potential (GRECP) [22| 
simulating interaction with the explicitly excluded Is to 
4/ electrons of Hf is used in 20-electron calculations of 
HfF + . In 10-electron calculations which are substantially 
less time consuming, 5s and bp spinors of hafnium and 
Is orbital of fluorine are frozen from the states averaged 
over the nonrelativistic configurations 5d 2 6s 6 6p 0A for 
Hf + and 2s 2 2p 5 for F and not treated explicitly. The 
generalized correlation atomic basis set (GCBS) 2 



24] 



(12sl6pl6dl0/10g)/[6s5p5d3/lg] is constructed for Hf. 
The ANO-L (14s9p4d3/)/[4s3p2dl£J atomic basis set 
listed in the MOLCAS 4.1 library [25j was used for flu- 
orine. The molecular orbitals are obtained by the com- 
plete active space self-consistent field (CASSCF) method 
[H, HI] with the spin-averaged part of the GRECP [13], 
i.e. only scalar-relativistic effects are taken into account 
at this stage. In the CASSCF method, orbitals are sub- 
divided into three groups: inactive, active and virtual. 
Inactive orbitals are doubly occupied in all the config- 
urations, all possible occupations are allowed for active 
orbitals, whereas virtual orbitals are not occupied. So, 
wave function is constructed as a full configuration inter- 
action expansion in the space of active orbitals and both 
active and inactive orbitals are optimized for subsequent 
correlation calculations of HfF + . Using the C 2v point 
group classification scheme, five orbitals in Ai, four in 
Bi and B2 and two in A2 irreps are included into the 
active space. In 10-electron calculations, one orbital in 
the Ai irreps (which is mainly 2s orbital of F) belongs to 
the inactive space. In 20-electron CASSCF calculations, 
the 5s and hp orbitals of Hf and Is orbital of F are added 
to the space of inactive orbitals. 

Next, the spin-orbit direct configuration interaction 
(SODCI) approach [H, El, HoJ modified by our group 
to account for spin-orbit interaction in configuration se- 
lection procedure [3l| with the selected single and dou- 
ble excitations from some multiconfigurational reference 
states ( "mains" ) is employed on the sets of different AS 
many-electron spin- and space-symmetry adapted basis 
functions (SAFs). Details on features of constructing 
the reference space and selection procedure are given in 
Refs. [H^il. 

Ten lowest states with the leading configura- 



tions 



( 3 Ai, 2 , 3 ; 1 



\..\a\o-\-K l ( 3 n -. +.i,2; 1 n) were calculated. Here 
<7i orbital is mainly formed by 2p z orbital of F with 
admixture of 6p z and 6s orbitals of Hf, 02 is mainly 6s 



orbital of Hf with admixture of 6p z orbital of Hf, 5 and 
7r are mainly the bd orbitals of Hf. 

To obtain spectroscopic parameters, six points between 
3.1 and 4.0 a.u. and point at 100 a.u. of the HfF + poten- 
tial curves were calculated for ten lowest-lying states in 
10-electron calculations and for four states in 20-electron 
ones. The 20-electron calculation is substantially more 
time-consuming and for the rest six states were calcu- 
lated only for one point, 3.4 a.u., in the present study. 
Comparing the latter calculations with corresponding 10- 
electron ones the core (5s 2 , bp 6 shells of Hf and Is 2 shell 
of F) relaxation and correlation corrections to the T e val- 
ues, called "core corrections" below and in Table HI were 
estimated. The core properties A\\ and -E e ff were calcu- 
lated only for the working 3 Ai state at point 3.4 a.u., 
which is close to the equilibrium distance (see Table JJ. 
Before calculating core properties the shapes of the four- 
component molecular spinors are restored in the inner 
core region after the two-component GRECP calculation 
of the molecule. For this purpose the nonvariational one- 
center restoration (NOCR) method 0, El 0, ED, EE El 
is applied. 

Results and discussion. The results of calculations for 
HfF + spectroscopic parameters are presented in Table HI 
The first point to note is that the cation is deeply bound. 
Second important result is that the 1 £ + state appears to 
be the ground one and the working state, 3 Ai, is the first 
excited one. Besides, the excitation energy from 1 E + to 
3 Ai is increased from 866 cmT 1 to 1633 cmT 1 after in- 
cluding the 5s and bp shells of Hf and Is shell of F to 
the relativistic correlation calculation. Excitation ener- 
gies from 1 £ + to other calculated low-lying states are 
also increased. Note also that the values obtained in 10- 
electron calculations for lowest four states when just the 
core correction described above is taken into account are 
in a good agreement with the purely 20-electron calcu- 
lations. (Note that accounting for correlation/relaxation 
of the 4/-shell would be also desirable but it is too con- 
suming and we expect it will result in smaller change of 
spectroscopic properties than that for the above core cor- 
rection.) The SO splittings of the 3 A and 3 n states are 
mainly due to the SO splitting of the 5c? shell of Hf: 



H 



(l5d'S5d) 



(3) 



The atomic Dirac-Fock calculation of 

a = 1269 cm" 1 



Hf + gives 

£5/2— £3/2 = 3173 cm -1 => a = 1269 cm -1 , where £ 5 / 2 
and £3/2 are orbital energies of hd§/2 and 5g?3/2 states of 
Hf + . The SO interaction ([3|) averaged over the 3 A or 3 n 
states is reduced to 



H 



LS 



A-(L-S), 



(4) 



where A = 1269/2 = 635 cm -1 , L and S are the orbital 
and spin momenta of HfF + . The SO interaction (j4|) leads 
to splitting between components of the 3 A and 3 n states 
on 1269 and 635 cm -1 , correspondingly. It is in a good 
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TABLE I: Calculated spectroscopic parameters for HfF 



State R e A 


T e cm' 1 


T e with core 


We cm 1 


D e cm 1 








correction" 






10-electron calculation 


X E+ 


1.784 








751 


51107 


3 Ar 


1.810 


866 


1599 


718 




3 A 2 


1.809 


1821 


2807 


719 




3 A 3 


1.807 


3201 


4324 


721 




X A 2 


1.814 


9246 


11519 


696 




X- 


1.856 


9466 


11910 


689 




3 n 0+ 


1.854 


9753 


12196 


699 




3 nr 


1.860 


10190 


12686 


687 




3 n 2 


1.856 


11898 


14438 


703 






1.870 


12642 


14784 


679 




20-electron calculation 


X E+ 


1.781 







790 


51685 


3 Ar 


1.806 


1633 




746 




3 A 2 


1.805 


2828 




748 




3 A 3 


1.804 


4273 




749 





a See paragraph "Methods and calculations" for details. 



agreement with the splitting of the 3 A state calculated 
in Table [T] but not with the 3 II one because of the SO 
interaction with closely lying 1 A and 1 IT states. 

The calculated Aii and E c g for the 3 Ai state are pre- 
sented in Table HH In opposite to A\\ [F], A\\ [Hf] and E e g 
are not seriously changed when the outer core electrons 
are included into calculation. Such behavior of Au\F] is 
explained by the fact that fluorine can be considered with 
good accuracy as a closed shell subsystem (thus having 
An [F] ps 0) in the 3 A state, i.e. large compensation of 
contributions from orbitals with different projections of 
total electronic momentum for An[F] takes place. There- 
fore, even small perturbation can seriously change the 
Ay [F] value. Calculated E c g is large and comparable 
with the corresponding value for the a(l) state of the 
PbO molecule [2l| . If one does not pay attention to the 
difference of signs our value is in 1.34 times larger than 
value obtained in the paper in scalar-relativistic cal- 
culations. Our rather crude estimate of the 3 Ai lifetime 
for the 3 Ai — > 1 £ + transition is about 1/2 sec. This value 
is difficult for accurate calculations because of small ab- 
solute values of both transition energies and, particularly, 
transition dipolc moments between those states whereas 
absolute errors are similar to those for other transitions. 
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TABLE II: Calculated parameters Ay [Hf] and A\\[F] (in 
MHz) and E cft (in 10 24 Hz/(e • cm)) for the 3 Ai state of 
177 Hf 19 F + at internuclear distance of 3.4 a.u. 



A,|[Hf] ^4 1| [F] E cB 
10-electron calculation 



-1250 -33.9 5.89 



20-electron calculation 



-1239 -58.1 5.84 
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